Increasing entanglement between Gaussian states by coherent photon subtraction 
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We experimentally demonstrate that the entanglement between Gaussian entangled states can 
be increased by non- Gaussian operations. Coherent subtraction of single photons from Gaussian 
quadrature-entangled light pulses, created by a nondegenerate parametric amplifier, produces de- 
localized states with negative Wigner functions and complex structures more entangled than the 
initial states in terms of negativity. The experimental results are in very good agreement with the 
theoretical predictions. 
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PACS numbers: : 03.67.-a, 03.65.Wj, 42.50.Dv 

Entanglement plays a key role in quantum information 
processing (QIP). Entanglement distillation [1], demon- 
strated for discrete- variable systems (ebits) in recent ex- 
periments 0, 0, Q| , allows one to produce strong entangle- 
ment between distant sites, initially sharing a larger set 
of weakly entangled states, and constitutes the basis of 
quantum repeaters, essential for long-distance quantum 
communications. An interesting alternative to discrete- 
level systems are quantum continuous variables (QCVs). 
In this case the information is encoded in the quadra- 
tures x and p of traveling light fields, which can be effi- 
ciently measured by homodyne detection. Optical para- 
metric amplification allows one to produce quadrature- 
entangled beams, used in many QIP protocols. To- 
gether with linear optics, these tools preserve the Gaus- 
sian character of the states involved in most of QCV ex- 
periments : the quasi-distributions (Wigner functions) 
of their quadratures remain Gaussian. However, it has 
been shown that Gaussian entanglement distillation re- 
quires non- Gaussian operations [H, @, 0]. Among sev- 
eral proposals [8|, |9|, one of the simplest is the condi- 
tional subtraction of photons from Gaussian entangled 
beams [Io|, 

EH [Hill, h y reflectin g 

a small part of these 
beams towards two photon-counting avalanche photodi- 
odes (APDs). If the reflectivity is low, a simultaneous 
detection of photons by the APDs heralds the subtrac- 
tion of exactly one photon from each beam. Recently, 
such methods allowed the preparation and analysis of 
several states with negative Wigner functions, including 
one- and two-photon Fock states , EliEfij], delocalized 



single photons [17|, LL8| and photon-subtracted squeezed 
states, very similar to quantum superpositions of coher- 
ent states with small amplitudes [13, [20[ . 

In this Letter, we experimentally demonstrate that 
non- Gaussian operations allow us to increase the entan- 
glement between Gaussian states, with a protocol pre- 
sented on Fig. [U An optical parametric amplifier (OPA) 
produces Gaussian quadrature-entangled light pulses, 
known as two- mode squeezed states [2l|. We pick off 
small fractions of these beams, which interfere with a 
well-defined phase on a 50/50 beam splitter (BS), and 
we detect photons in one of the BS outputs. This way, 



we subtract a single photon delocalized in the two beams 
and prepare a complex quantum state with a negative 
two-mode Wigner function. We determine a range of 
experimental parameters where the entanglement of the 
prepared state, quantified by the negativity [22[, is sig- 
nificantly higher compared to that of the initial Gaussian 
state. 

Operating with single photon counts rather than with 

coincidences as proposed e. g. in m m m mm, 

this protocol allows for much higher generation rates and 
produces states more robust to experimental imperfec- 
tions (see below). Besides, it is more efficient at moderate 
OPA gain: in the zero-gain limit, the detection of a pho- 
ton transforms a state with almost no entanglement into 
a maximally entangled ebit state ( 1 10) + \01))/y/2. With 
the higher gain (up to 3 dB) used in the present experi- 
ment, the generated states have a much richer structure 
as is shown below. 

Our experimental setup is presented in Fig. O Nearly 
Fourier-limited femtosecond pulses (180 fs, 40 nJ), pro- 
duced by a Tiisapphire laser with a 800 kHz repeti- 
tion rate, are frequency doubled by a single pass in a 
100 jam-thick type I noncritically phase-matched potas- 
sium niobate (KNbOa) crystal. The frequency-doubled 
beam pumps an identical crystal used as an optical para- 
metric amplifier (OPA), generating Gaussian quadrature- 
entangled pulses spatially separated by an angle of 10°. 
Adjusting the pump power allows us to vary the two- 
mode squeezing between and 3.5 dB. The photon pickoff 
beam splitters are realized with a single polarizing beam 
splitter (PBS) cube, where the signal and idler beams 
are recombined spatially but remain separated in polar- 
ization. A small adjustable fraction R of both beams 
is sent into the APD channel, where they interfere on a 
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FIG. 1: Coherent photon subtraction from Gaussian entan- 
gled beams. 
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FIG. 2: Experimental setup. 



50/50 BS. A tilted half- wave plate compensates for resid- 
ual birefringence. An APD detects one of the 50/50 BS 
outputs after spatial and spectral filtering. The signal 
and idler beams transmitted through the pickoff beam 
splitter are projected into a non-Gaussian state by an 
APD detection. They are spatially separated on another 
PBS, where they are combined with bright local oscilla- 
tor beams. A quarter- wave and a half- wave plate allow 
us to prepare two local oscillators with equal intensities 
and a well-defined relative phase. The signal and idler 
beams are analyzed by two time-resolved homodyne de- 
tections, which sample each individual pulse, measuring 
one quadrature £1,2 (#1,2) m phase with the local oscilla- 
tor. 

In this setup, all the relative phases except <pi and 
<p2 (see Fig. [2j) are precisely adjusted with wave plates. 
Phase fluctuations concern only the initial two-mode 
squeezed state, where the phase difference is not defined 
and plays no role, and the slow (thermal and acous- 
tic) phase sum fluctuations simply rotate the two-mode 
squeezing ellipse and can be compensated by shifting the 
common phase fa of the local oscillators. This phase can 
be scanned with a piezo translator, and rapidly measured 
using the unconditioned two- mode squeezing variance. 

Quantum states with negative Wigner functions are 
very sensitive to experimental imperfections. In our 
case, the most important issue are spurious APD trig- 
ger events, due to imperfect filtering, limited qualities of 
the optical beams, imperfect mode-matching between the 
subtracted beams, and APD dark counts. An APD count 
corresponds to the desired subtraction event with a suc- 
cess probability £ < 1. This explains why single-photon 
protocols are more robust than two-photon ones, where 
the total success probability is only £ 2 . Another issue 
is the OPA excess noise. To describe it, we can consider 
that a first amplification process creates a pure entangled 



state with a two-mode squeezing variance s 



, and 



that each of the resulting modes is independently ampli- 
fied with a gain h = cosher) by a phase-independent 
amplifier with a relative efficiency 7. The finite homo- 
dyne efficiency rj and the homodyne excess noise e also 



deteriorate the measured data. However, they are not in- 
volved in the generation process but only in the analysis, 
and we can correct for their effects in order to determine 
the actual Wigner function of the generated state. Even 
with none of these imperfections, this protocol would still 
be limited by the finite pick-off BS reflectivity i?, required 
for a sufficient APD count rate but inducing losses on the 
transmitted beam. The limited overall efficiency ji — 5% 
of the APD channel has little effect in this experiment. 

A detailed analytic model including all these 

imperfections yields an expression for the Wigner func- 
tion W of the state studied in our experiment : 

W(x u p 1: x 2 ,p 2 ) = W 3 (x+,p+)W c (x-,p-) (1) 



where x± = Xl ^ 2 , P± = Pl ^ 2 ? W a is the Wigner func- 
tion of a single-mode squeezed state, and W c corresponds 
to a photon-subtracted squeezed state analyzed in [l9[. 
More explicitly : 
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In this experiment the photon-subtracted state is "delo- 
calized" into two spatially separated modes 1 and 2 and 
revealed by measuring the correlations between identical 
quadratures, the ant icorrelat ions remaining in the initial 
squeezed state. 

Without assuming any particular shape for W s and 
W c , we can experimentally show that the state becomes 
separable if we make a joint measurement, transform- 
ing xi ? 2 into x± by rotating the polarizations by 45° 
with the optional half- wave plate shown on Fig. O We 
then observe that the quadratures measured by one de- 
tection do not depend on the other (see Fig. [3]). For 
every 0± the joint distribution, and hence the Wigner 
function, becomes factorable : P(a; + (# + ), X- (0-)) = 
P s (x+(0+))P c (x-(0-)). It means that one can fix 0- = 
6+ = and scan to perform a complete tomography of 
this state. 

This property considerably simplifies the experimental 
analysis. We perform direct homodyne measurements of 
the entangled quadratures xi^(0), keeping the entangled 
modes 1 and 2 separated without mixing them. We use 
the fact that the state factorizes in the x + , X- basis to re- 
construct it from a limited set of data : instead of a time- 
consuming general two-mode tomography, which requires 
to measure x\{6i), ^2(^2) with all possible combinations 
of phases, we can restrict ourselves to 9\ = 2 = 6. In 
practice, we set the relative phase between the local os- 
cillators to zero, scan the common phase ^2 = 0, measure 
#1,2(0) an d calculate We reconstruct the distri- 

butions P c (x-(0)) and P 3 (x+(9)) for several phases. We 



FIG. 3: State separability test after interference between sig- 
nal and idler beams : (a) Joint distribution P(x+,x_), (b) 
Distributions P s (x+) and P c (x_), for 11 values of X- (resp. 

chosen between -2 and +2. This separability was ver- 
ified for several randomly chosen values of 6+ and 6- (here 
0+ = 20° and 0- =50°). 



observe that the measured distributions are invariant un- 
der 6 — > 7r±#, so we restrict the analysis to < 6 < tt/2. 
Typically, we measure 6 to 12 different quadrature dis- 
tributions, with 10000 to 20000 data points each. A nu- 
meric Radon transform allows us to reconstruct the un- 
corrected Wigner functions W c and W s . We can correct 
for the homodyne detection losses (77 = 70%, e = 1% 
of the shot noise) using a maximal-likelihood algorithm 
25l . [2^ to obtain the Wigner function W of the generated 
state. We use W to calculate the density matrix p of the 
state and obtain its entanglement, given by the negativ- 
IIp Ti I 



ity Af 



— , where T\ is the partial transposition 



operation [22[. 

Figure [4] presents the tomography of a state produced 
with 1.8 dB of squeezing and a BS reflectivity R = 5%. 
The Wigner function, corrected for detection losses, is 
clearly negative : W c (0) = -0.13 ± 0.01 (0.01 ± 0.01 
before correction). The entanglement of this state is Af = 
0.34 ±0.02, whereas for the initial state (before the pickoff 
BS) M) = 0.24 ±0.01. 

In Refs. 

he mi 

we demonstrated another analysis 
method, more constrained but also much faster and closer 
to the physics of the experiment. If we assume that the 
Wigner function has the form defined in Eq. [TJ we can 
easily extract the parameters a, A, 6, B from the second 
and fourth moments of the measured distributions, and 
determine the Wigner function, the density matrix, and 
the quadrature distributions of the measured state. For 
a given squeezing, one can also obtain from a, A, b and 
B the values of all the experimental parameters intro- 
duced above. To correct for homodyne losses, we sim- 
ply calculate the Wigner function that we would mea- 
sure with an ideal detection (77 = 1 and e = 0), using 
the values extracted from the experimental data for all 
the other parameters. As shown in Fig. [H the distri- 
butions reconstructed with this method are in excellent 
agreement with those directly extracted from the data, 
and the Wigner function is almost indistinguishable from 
the one obtained with the maximal-likelihood algorithm. 
Both methods give the same values for the negativity. 
A natural question to ask is whether this protocol 
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FIG. 4: (a) Set of experimentally measured quadrature dis- 
tributions (dots), compared to those reconstructed from our 
model (solid line), (b) Wigner function corrected for homo- 
dyne detection losses, obtained with a standard maximal- 
mikelihood algorithm (MaxLike), compared to the result of 
our model. This state is produced with 1.8 dB of squeezing 
and R 5%. 



works for an arbitrary squeezing or if, when the initial 
state is already strongly entangled, by performing an im- 
perfect photon subtraction we actually lose more entan- 
glement than we gain. It has already been shown that, 
when the pick-off beam splitters have a finite reflectivity, 
subtracting one photon from each of the Gaussian en- 
tangled beams may actually decrease the entanglement 
[10[. Using our model, we can take into account all the 
other experimental parameters to derive an analytic ex- 
pression for p Tl , which can be diagonalized numerically 
in a few seconds to obtain the expected negativity of a 
given state. We found that the experimental imperfec- 
tions have a very strong effect. For example, for an initial 
squeezing of 3 dB, the negativity increases ideally from 
M) 0.50 to Af = 0.90. If we assume R 3% for the 
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FIG. 5: Tomography of two states produced with 1.3 dB (first 
row) and 3.2 dB (second row) of squeezing and R = 10%: 
three different cuts of the two-mode Wigner function. 
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FIG. 6: Entanglement negativity of the initial and final states 
as a function of squeezing for several pickoff BS reflectivities, 
corrected for homodyne detection losses. Solid lines are theo- 
retical calculations using the average values of the experimen- 
tal parameters. 

pick-off BS, then M = 0.81, but if we include the average 
values of experimental parameters involved in the state 
preparation, (7) = 0.22 and (£) = 0.78, M drops down to 
0.51, whereas for the initial state, only slightly affected 
by 7, Afo = 0.49. 

To verify this experimentally, we performed several 
tomographies for different BS reflectivities and degrees 
of squeezing. Figure [5] presents two tomographies with 
R = 10% for a small (1.3 dB) and a high (3.2 dB) squeez- 
ing. As expected, the Wigner function becomes more 
phase-dependent and less negative as the state becomes 
"bigger" and more sensitive to decoherence. 

Figure [6] shows the entanglement negativity of the 
photon-subtracted states and the corresponding initial 
states. The solid lines are theoretical calculations us- 
ing the average values of the experimental parameters 
involved in the state preparation, (7) =0.22 and (£) = 
0.78. Two domains appear on the graph: the upper left, 
where this process actually increases the entanglement, 
and the lower right, where the initial state is too sen- 
sitive to the added losses. As expected, this protocol is 
particularly efficient at low squeezing. One can show that 
when the squeezing and hence the entanglement of the 
initial state tend to 0, the negativity of the final state 

has a nonzero limit A/" r ^o = ^° p - °^ where 
C = ^iq~ir • At low squeezing, experimental imperfec- 
tions have a moderate effect on the state. We never- 
theless succeed in improving the negativity of a state 
with up to 3 dB of squeezing, which corresponds to a 
strong squeezing regime where small experimental im- 
provements strongly affect the performance of the pro- 
tocol. For example, increasing £ by a mere 4% with 
R = 3% should displace the crossover point from 3 to 
4 dB (in other experiments, where mode matching be- 



tween subtracted photons was not an issue, £ reached 0.9 

0). 

In conclusion, the present photon subtraction protocol 
allows one to increase the entanglement between Gaus- 
sian states with up to 3 dB of squeezing, and even small 
experimental improvements should significantly increase 
this limit. For QIP protocols specifically requiring Gaus- 
sian entanglement, these non-Gaussian states could in 
principle be used as a starting point for a "Gaussifica- 
tion" procedure @. This demonstrates one of the key 
steps required for long distance quantum communica- 
tions with continuous variables. 
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